For  Reference 


NOT  TO  BE  TAKEN  FROM  THIS  ROOM 


Six  UBBtf 
attmamis 

t 


University  of  /Wbet'fn 
Printing  Department 


Digitized  by  the  Internet  Archive 
in  2018  with  funding  from 
University  of  Alberta  Libraries 


https://archive.org/details/effectofmagneticOOralp 


UNIVERSITY  OF  ALBERTA 


in  4 


THE  EFFECT  OF  A  MAGNETIC  FIELD 

ON 

NORMAL  SEGREGATION  IN  ALUMINUM(P ICH) -  COPPER  ALLOYS 

by 

Ralph  Clarence  Dorward 


A  THESIS 

SUBMITTED  TO  THE  FACULTY  OF  GRADUATE  STUDIES 
IN  PARTIAL  FULFILMENT  OF  THE  REQUIREMENTS 
FOR  THE  DEGREE  OF  MASTER  OF  SCIENCE 

DEPARTMENT  OF  MINING  AND  METALLURGY 

EDMONTON,  ALBERTA 


APRIL  1964 


ABSTRACT 


The  effect  of  a  magnetic  field  on  the  solidification  of  aluminum 
(rich)  -  copper  alloys  was  investigated  using  normal  segregation  as  a 
criterion.  The  field  strength  was  34,000  oersteds  applied  perpendicular 
to  the  solidification  direction  for  freezing  rates  ranging  from  1/4  to  2 
inches  per  hour.  For  a  solidification  velocity  of  l  to  2  inches  per  hour 
there  was  an  increase  of  12  pet.  in  the  segregation  for  the  "field"  solid¬ 
ified  ingots  compared  to  that  for  the  normally  solidified  ingots.  However, 
for  freezing  rates  of  1/2  to  1/4  inches  per  hour  the  field  did  not  affect 
the  segregation  to  any  observable  extent.  A  theory  has  been  presented 
which  qualitatively  explains  the  results. 
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I.  INTRODUCTION 

The  electronic  structure  and  transport  properties  of  solids, 

such  as  galvanomagnetism,  thermomagnetism,  and  the  Hall  effect  have 

been  successfully  studied  by  the  application  of  a  magnetic  field.  This  leads 

to  the  hypothesis  that  a  magnetic  field  may  affect  solidification  phenomena 

either  through  its  effect  on  the  equilibrium  phase  relationship  of  the  alloy 

system, or  through  the  irreversible  transport  processes  of  solidification. 

A  study  of  solidification  phenomena  in  a  strong  magnetic  field  should, 

therefore,  elucidate  crystallization  mechanisms  and  the  nature  (electronic) 

of  the  solid  -  liquid  system,  provided  the  theoretical  principles  of  alloy 

solidification  and  magnetic  field  interactions  with  electrically  conducting 

fluids  (magnetohydrodynamics)  are  correctly  applied. 

Previous  work  by  D.  R.  Colton^  ^  on  the  solidification  of  binary 

alloys  of  bismuth  with  antimony  and  tin  in  a  13,000  oersted  magnetic  field, 

applied  parallel  to  the  direction  of  solidification,  indicated  that  the  field 

changed  the  amount  of  inverse  segregation  by  about  4*100  pet.  and  —70  pet. 

(2) 

respectively  over  the  "no  -  field"  results.  J.  R.  Cahoon  investigated 
the  effect  of  a  magnetic  field  on  inverse  segregation  in  the  aluminum  (rich)  - 

copper  system,  and  found  an  increase  in  segregation  of  60  pet.  for  a 

/ 

magnetic  field  strength  of  34,000  oersteds.  The  magnitude  of  segregation 
was  found  to  be  independent  of  field  direction.  The  observed  segregation 
in  the  "field"  solidified  ingots  for  each  system  can  be  theoretically 
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calculated  by  the  method  of  Kirkaldy  and  Youdelis  if  the  equilibrium 
phase  diagram  is  altered  by  an  appropriate  shift  in  the  liquidus  curve, 
which  is  due  to  a  proposed  diffusion  inhibition  in  the  liquid  caused  by  the 
field.  In  order  to  test  further  the  validity  of  the  alteration  of  the  cons¬ 
titutional  diagram  it  was  decided  to  study  the  effect  of  a  magnetic  field  on 
normal  segregation  in  the  aluminum  (rich)  -  copper  system,  since  this 
would  allow  better  control  (and  evaluation)  of  the  solidification  variables. 
The  amount  of  segregation,  which  is  readily  determined  by  ordinary 
chemical  analysis,  is  related  to  the  equilibrium  phase  relationship  of  the 
system.  Any  change  in  this  relationship  would  be  revealed  by  a  change  in 
the  normal  segregation  pattern  of  a  directionally  solidified  ingot. 
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II.  HISTORICAL  REVIEW 

A.  Theoretical  Considerations 
Normal  freezing  is  a  process  of  controlled  solidification  in 
which  the  entire  body  under  consideration  is  melted  and  then  frozen  slowly 
and  progressively  from  one  end.  Complications  such  as  dendrite  formation 
and  inverse  segregation  present  under  conditions  of  rapid  solidification 
(casting)  are  absent.  In  the  following  a  brief  resume  of  the  theory  of  normal 
solidification  is  given. 

(1)  The  Distribution  Coefficient 

Considering  the  partial  equilibrium  diagram  shown  in  Figure  1, 
solidification  of  a  binary  liquid  alloy  of  solute  concentration  C0  begins  when 
the  temperature  reaches  T0.  The  concentration  of  the  first  solid  to  form 
is  given  by  the  intersection  of  the  horizontal  at  temperature  T0  with  the 
solidus.  It  is  convenient  to  introduce  an  equilibrium  distribution  coefficient, 
Kc,  which  is  defined  as  the  ratio  of  the  equilibrium  concentration  of  the 
solute  in  the  solid  (at  the  solid  -  liquid  interface)  to  the  equilibrium  con¬ 
centration  of  the  solute  in  the  liquid.  The  distribution  coefficient  may  be 
greater  or  less  than  unity  depending  on  whether  the  solute  raises  or  lowers 
the  melting  point  of  the  solvent.  Figure  1  illustrates  the  case  for  K0less 
than  unity;  whereas  if  the  liquidus  had  a  positive  slope,  K0  would  be  greater 
than  unity.  The  distribution  coefficient  may  be  related  in  part  to  the  dis¬ 
tortion  the  solute  atoms  cause  in  the  solid  lattice.  For  Kc  equal  to  unity 


-  -J  .) 


. 


Figure  1.  Partial  Equilibrium  Diagram  with  K<1 


Figure  2.  Concentration  Gradient  in  the  Liquid  Ahead  of  an 


Advancing  Solid  -  Liquid  Interface 


the  solute  creates  no  more  distortion  in  the  solid  lattice  than  it  does  in  the 
liquid,  whereas  if  it  is  less  than  unity  the  solute  lowers  the  melting  temp¬ 
erature  of  the  solid,  thereby  weakening  the  bonding  between  the  atoms  and 
increasing  the  interatomic  distance.  The  opposite  is  true  for  KQ  greater 
than  unity.  The  departure  of  the  distribution  coefficient  from  unity  is  a 
measure  of  the  degree  of  lattice  distortion  in  a  binary  alloy  system,  assuming 
no  electronic  (valency  or  electro -chemical)  effects  are  present. 

Macroscopically  the  equilibrium  distribution  coefficient  is  rarely 
achieved  in  a  real  freezing  process,  for  the  reason  that  equilibrium  con¬ 
ditions  are  not  present  (macroscopically)  even  for  very  slow  solidification 

rates.  Microscopically,  however,  local  equilibrium  is  assumed  to  exist 

(3) 

at  the  interface.  Consider  the  case  for  Kcless  than  unity.  This  corres¬ 
ponds  to  "rejection"  of  solute  atoms  by  the  solid  as  it  forms,  and  therefore 
the  solute  concentration  in  the  remaining  liquid  increases.  This  situation 
is  illustrated  in  Figure  2,  which  shows  the  concentration  profile  near  an 
advancing  solid  -  liquid  interface.  At  the  interface  the  concentration  freezing 
out  in  the  solid  is  Ko  times  the  concentration  in  the  liquid.  If  diffusion  is  the 
only  transport  mechanism  in  the  liquid  there  are  two  factors  which  influence 
the  forms  of  the  concentration  gradient  ahead  of  the  interface: 

(  i)  diffusion  in  the  liquid  tending  to  spread  the  solute  throughout  the  melt 
at  a  rate  determined  by  the  diffusion  coefficient  D,  and 
(ii)  the  rate  of  freezing  R,  which  acts  as  a  source  of  solute  at  the  interface. 
When  the  solute  concentration  in  the  interfacial  liquid  reaches  the  value 
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C0/K0,  steady  state  conditions  exist  because  the  concentration  Cc  of  the 

liquid  entering  the  "solidification  zone"  equals  that  leaving  it  in  the  solid. 

(4) 

Tiller  et  al  derive  the  following  expression  for  the  concentration  gradient 
for  the  steady  state  condition: 


^  u 


H- 


l  -  Ko 

K0 


(1) 


where  x'  is  the  distance  from  the  interface  into  the  liquid. 

Diffusion,  however,  is  not  the  only  transport  mechanism  in  the 
liquid.  There  is  usually  some  degree  of  convective  mixing  even  in  unstirred 
melts  and  a  steady  state  condition  is  not  attained.  Any  turbulence  in  the 
liquid  affects  the  concentration  gradient  ahead  of  the  solid  -  liquid  interface 
by  reducing  its  thickness.  Thus,  the  solute  concentration  in  the  liquid  at 
the  interface  will  be  lowered  from  that  shown  in  Figure  2;  and  if  solid¬ 
ification  rates  are  sufficiently  slow,  the  liquid  will  attain  a  uniform  concen¬ 
tration  through  diffusion  and  convective  mixing. 


(2)  Normal  Segregation 

As  solidification  proceeds  along  the  charge,  rejected  solute 
accumulates  in  the  liquid.  Thus,  due  to  this  effect,  and  also  because  the 
liquid  volume  is  being  decreased,  the  liquid  solute  concentration  rises. 
This,  in  turn,  causes  the  concentration  freezing  out  to  increase  since  it 
is  K0  times  that  in  the  liquid.  If  freezing  is  considered  to  proceed  along 
a  cylinder  from  one  end  to  the  other,  an  end  -  to  -  end  segregation  will 
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be  produced  in  the  solid  ingot.  This  is  known  as  normal  segregation. 

Depending  on  the  conditions  of  freezing  four  different  segregation 
forms  can  result  as  illustrated  in  Figure  3.  They  are: 

(i)  If  solidification  occurs  slowly  enough  for  equilibrium  to  be  maintained 

there  will  be  no  segregation  in  the  solid  because  time  is  allowed  for  complete 

-"I 

diffusion.  Diffusion  rates  in  solid  metals  are  so  slow  (D  '*"10  cm.  /sec.), 
however,  that  equilibrium  freezing  almost  never  occurs  in  practice. 

(ii)  When  solidification  proceeds  slowly  enough  for  diffusion  in  the  liquid 

-5  2. 

(D  10  cm.  /sec.)  to  be  complete  but  fast  enough  for  diffusion  in  the  solid 

(5) 

to  be  negligible,  Pfann  shows  that  the  concentration  gradient  in  the  frozen 
solid  is  given  by 

C5  =  K0C0  (i  -%)  K°  ,  (2) 


where  is  the  fraction  solidified.  Under  these  conditions  the  maximum 
possible  overall  segregation  is  obtained. 

(iii)  When  the  freezing  rate  is  rapid  compared  with  diffusion  and  convective 

processes  in  the  liquid,  the  distribution  coefficient  rises  to  unity  and  a 

plateau  is  produced  in  the  segregation  curve  (cf.  Equation  2).  Under 

,  (4) 

conditions  of  negligible  convection,  Tiller  et  al  show  that  the  solute 
concentration  in  the  solid  for  the  initial  transient  is  given  by 


CS  =  Co  (\-  ke) 


-  e 


“  D  X 


+■  K, 


(3) 


where  x  is  the  distance  measured  from  the  starting  point  of  solidification. 
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Solute  Concentration 
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(i)  Equilibrium  Freezing 

(ii)  Complete  diffusion  in  the  liquid  and  negligible  diffusion 
in  the  solid 

(iii)  Convective  mixing  in  the  liquid  insignificant 

(iv)  Interfacial  diffusion  layer  present  with  convective  mixing 
in  the  liquid 

Figure  3.  Segregation  Patterns  Corresponding  to  Various 
Solidification  Conditions. 
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This  type  of  normal  segregation  is  not  often  encountered  in  actual  solid¬ 
ification  because  the  freezing  rates  required  are  difficult  to  maintain 
without  resulting  in  dendritic  growth  and  inverse  segregation. 

(iv)  In  general  solidification  conditions  are  somewhere  between  (ii)  and 
(iii),  ie.  the  liquid  is  neither  homogenous  as  assumed  in  Pfann's  Equation 
(2),  nor  cored  to  the  extend  expressed  in  Equation  (l).  This  is  because 
convective  mixing  will  occur  in  the  liquid,  and  it  will  be  homogenous  except 
for  a  stagnant  surface  layer  of  average  thickness  <$.  Burton  et  al  ^  have 
shown  that  for  this  case  Pfann's  equation  can  still  be  used  for  the  concen¬ 
tration  distribution  in  the  solid,  but  now  an  effective  distribution  coefficient, 
K,  must  be  used  in  place  of  K0  to  take  account  of  the  effect  of  the  stagnant 
unmixed  layer.  They  derive  the  effective  distribution  coefficient  as: 

K  =  - -  K\  •  <4> 

k0  +  0-  K„)  e  0 

Their  derivation  requires  a  planar  interface,  and  complete  mixing  in  the 

liquid  alloy  is  assumed  except  for  the  interfacial  layer  of  constant  average 

thickness.  The  thickness  of  this  layer  depends  on  the  efficiency  of  mixing 

in  the  bulk  liquid,  and  transport  of  solute  through  it  takes  place  only  by 

(7) 

diffusion.  Dean  et  al  have  shown  that  when  the  rate  of  solidification  is 

-5 

relatively  slow  (in  the  order  of  10  cm.  per  second),  that  the  effective 
distribution  coefficient  approaches  the  equilibrium  value. 
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B,  Experimental  Observations 

During  the  past  few  years  a  great  deal  of  experimental  work 

has  been  conducted  in  the  field  of  normal  solidification.  In  the  following 

a  summary  is  given  of  some  of  the  work  substantiating  the  above  theories. 

The  existence  of  a  solute  concentration  gradient  ahead  of  an 

( 8 ) 

advancing  solid  -  liquid  interface  was  shown  by  Burton  et  al.  They 
determined  experimentally  the  effective  distribution  coefficients  of  various 
solute  elements  in  germanium  crystals  grown  from  the  melt.  By  varying 
the  thickness  of  the  layer  through  rotation  of  the  crystal,  they  found  that 
the  effective  distribution  coefficient  depended  on  growth  conditions  as 
expressed  in  Equation  (4). 

Using  the  expression  for  the  steady  state  concentration  gradient 
(4) 

(Equation  l),  Tiller  et  al  predicted  the  breakdown  of  a  macroscopically 

smooth  interface  at  a  critical  ratio  of  the  temperature  gradient  in  the 

(9) 

liquid  at  the  interface  to  the  rate  of  growth.  Walton  et  al  and  Tiller 

and  Rutter  experimentally  verified  this  theoretical  prediction  by 

examining  the  solid  -  liquid  interfaces  of  ingots  which  were  solidified  under 

controlled  conditions.  Convective  mixing  was  assumed  minimal  because 

Wagner  had  shown  that  this  effect  is  not  operable  in  the  systems  studied 

for  the  rates  of  solidification  used  (greater  than  0.  5  cm.  per  minute). 

(12) 

Y ue  and  Clark'  7  conducted  directional  freezing  experiments 
on  dilute  magnesium  alloys.  They  determined  the  effective  distribution 


coefficients  for  the  alloys  at  various  freezing  rates  and  discovered  that  the 
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normal  segregation  carves  followed  the  theoretical  expression  given  by 

Pfann  (Equation  2),  thus  showing  that  for  the  case  of  partial  mixing,  K6 

can  be  replaced  by  K.  They  also  show  that  at  very  slow  freezing  rates 

the  effective  distribution  coefficient  approaches  the  equilibrium  value. 

(7) 

Similar  experiments  by  Dean  et  al  have  also  led  to  the  same  general 
conclusions . 

(13) 

Weinberg'  determined  the  solute  distribution  in  the  liquid 
during  solidification  of  tin  -  silver  alloys  by  quenching  the  residual  liquid 
after  freezing  had  commenced  and  then  analyzing  the  solute  concentration 
ahead  of  the  "frozen  -  in"  interface.  He  found  that  when  convective  mixing 
in  the  liquid  had  a  negligible  effect, the  expressions  of  Tiller  et  al  (Equations 
l  and  3)  were  in  good  agreement  with  experiment.  He  also  found  that  when 
convective  mixing  affects  the  solidification  process,  the  segregation  in  the 
solid  may  be  expressed  by  Pfann's  equation,  providing  an  effective  dist¬ 


ribution  coefficient  is  used. 


(14) 


Recently  Niessen  et  al '  7  successfully  applied  Pfann’s  equation 


to  normally  solidified  ingots  of  zinc  -  gold  to  determine  the  effective 
distribution  coefficients  corresponding  to  various  growth  velocities.  They 
then  extrapolated  the  experimental  coefficients  to  zero  rate  of  growth  to 


obtain  the  equilibrium  distribution  coefficient  of  the  system. 
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III.  EXPERIMENTAL 

A.  Apparatus 


(l)  The  Magnet 

The  magnetic  field  for  the  solidification  experiments  was 

obtained  using  a  water  -  cooled)  pole  type  electromagnet  (see  Figure  4) 

(15) 

which  was  an  adaptation  of  the  Uppsala  Electromagnet'  .  The  magnet 
supplied  a  maximum  field  of  35,  000  oersteds  uniform  to  -  3  pet,  over  a 
cylindrical  volume  3.  5  inches  in  diameter  and  l.  5  inches  in  length.  A 
thirty  six  kilowatt,  direct  current,  constant  voltage  generator  which  gave 
a  maximum  voltage  of  48  volts  with  a  ripple  of  £2  pet.  supplied  power  to 
the  magnet. 


(2)  Solidification  Unit 

The  alloys  were  solidified  in  the  horizontal  stainless  steel 
tube  and  resistance  furnace  assembly  shown  in  the  photograph  (Figure  5). 

A  schematic  drawing  of  the  cross  section  of  the  unit  is  shown  in  Figure  7. 
The  tube  was  3  feet  long  and  3/4  of  an  inch  in  diameter,  and  was  suitably 
connected  to  provide  an  atmosphere  of  argon.  A.  copper  cooling  coil  was 
positioned  at  the  front  end  (where  solidification  was  initiated)  of  the  furnace 
in  order  to  provide  a  reproducible  temperature  gradient  at  the  solid  - 
liquid  interface.  The  alloy  ingots  were  contained  in  a  4  inch  long  cylind¬ 
rical,  recrystallized  alumina  boat  which  was  pushed  out  of  the  furnace  by 
a  moveable  stainless  steel  rod  1/4  of  an  inch  in  diameter.  The  rod  was 


- 


>v  '>  Jr  V  ,  '•  '  •  .  ;•>.  !.<•  .  ,  Oj  ' 


13 


Figure  4.  Photograph  of  the  Pole  -  type  Electromagnet. 


Figure  5.  Photograph  of  the  Solidification  Unit. 


Figure  6.  Photograph  of  the  Drive  Mechanism. 
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igure  7.  Cross  Section  of  the  Solidification  Unit 
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sealed  at  the  end  of  the  tube  with  an  "O"  ring.  Another  sealed  stainless 
steel  rod  of  the  same  dimensions  was  connected  to  the  front  end  of  the 
boat  for  mixing  purposes.  Steady  travel  rates  were  achieved  with  a 
synchronous  electric  motor  and  a  Zero  -  Max  (Figure  6)  connected  by  a 
flexible  drive  shaft  to  a  screw  mechanism.  The  temperature  was  controlled 
by  a  stationary  Chromel-P  -  Alumel  thermocouple  placed  at  the  rear  of 
the  furnace.  A  Wheelco  controller  and  a  high  -  low  variable  transformer 

,  O 

maintained  the  temperature  to  -  2  C. 

Initial  experiments  with  a  graphite  boat  indicated  that  the  alloy 

supercooled  when  withdrawn  from  the  furnace.  This  caused  the  solid  - 

liquid  interface  to  move  rapidly  along  the  first  inch  or  less  of  the  charge. 

Random  nucleation  occured  in  the  initial  portion  of  the  melt  resulting  in 

entrapment  of  the  solute  during  the  ensuing  rapid  freezing.  This  effect 

(7)  (12) 

has  also  been  reported  by  Dean  et  al  '  and  by  Y ue  and  Clark.  '  To 
largely  eliminate  random  nucleation  it  was  necessary  to  have  good  heat 
withdrawal  through  the  front  of  the  boat  and  low  heat  removal  through 
the  sides.  Recrystallized  alumina  tubes  with  graphite  end  plugs  were 
therefore  used  because  of  the  relatively  high  thermal  conductivity  of  graphite 
compared  to  alumina.  The  front  end  plug  was  also  tapered  in  the  manner 
described  by  Lawson  and  Nielsen'  ’  for  growing  single  crystals.  A.  copper 
block  was  fixed  to  the  leading  end  for  effective  heat  removal  where  solid¬ 


ification  was  initiated. 


>  I 

' 
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B.  Method 

The  alloys  were  prepared  in  alundum  coated  silica  crucibles 
by  melting  in  an  electric  furnace.  About  four  hundred  grams  of  alloy  con- 
taining  3  pet.  copper  were  prepared  in  each  melt  from  aluminum  (99.  91 
pet.  purity)  and  electrolytic  copper  (99.99  pet.  purity).  The  melts  were 
thoroughly  mixed  by  pouring  them  back  and  forth  between  two  crucibles 
ten  times.  After  mixing,  the  melts  were  degassed  by  bubbling  with  argon 
for  approximately  one  minute.  The  alloys  were  then  poured  at  a  superheat 
of  50°  C  into  a  cool  cylindrical  graphite  mould  with  dimensions  1  1/2  inches 
outside  diameter  and  5/8  inches  inside  diameter  by  8  inches  high.  Four 
ingots  were  obtained  from  each  melt.  A  cool  mould  was  used  so  that  solid¬ 
ification  would  proceed  quickly  from  the  sides  along  the  length  and  not  from 
the  top  or  bottom  in  the  manner  of  ingots  quickly  solidified  unidirectionally. 
However,  to  eliminate  any  possible  effects  of  inverse  segregation,  two 
inches  were  cropped  from  each  end  of  the  ingot  where  segregation,  if 
present,  would  be  greatest.  The  ingots  were  then  machined  down  to  fit 
into  the  alumina  boat. 

Before  an  experimental  run  the  boat  containing  the  alloy  rod 
was  positioned  in  the  stainless  steel  tube  so  that  the  solid  -  liquid  interface 
would  be  between  the  pole  pieces  of  the  magnet,  and  the  system  was  purged 
with  argon  for  twenty  minutes.  After  the  furnace  was  turned  on  and  the 


''  All  concentrations  are  in  weight  per  cent. 
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alloy  melted,  the  argon  atmosphere  was  kept  at  a  slight  positive  pressure 
for  the  duration  of  the  run.  Solidification  did  not  commence  until  time 
was  allowed  for  steady  -  state  temperature  gradients  to  become  established. 
During  this  time  the  boat  was  rotated  back  and  forth  at  five  minute  intervals 
by  means  of  the  stainless  steel  mixing  rod.  This  process  served  to  eliminate 
any  gravity  segregation  formed  on  melting,  and  furthermore,  it  assured  a 
uniform  solute  concentration  lengthwise.  The  alloy  was  then  solidified 
by  pushing  the  boat  through  the  furnace  at  the  desired  speed.  Care  was 
taken  to  reproduce  exactly  all  conditions  for  the  duplicate  runs,  one  in  a 
magnetic  field  of  34,000  oersteds  perpendicular  to  the  solidification  direc¬ 
tion  and  the  other  without  a  field. 

The  solidified  ingots  were  sectioned  and  the  swarfs  of  nine 
cuts  along  their  lengths  were  collected  for  analysis.  The  sections  were 

(17) 

retained  for  metallographic  examination.  The  iodide  method  of  analysis 
was  used  and  all  precautions  for  good  volumetric  procedure  were  followed 
to  ensure  good  precision  in  the  analysis.  Corresponding  slices  of  the 
duplicate  ingots  were  analyzed  at  the  same  time  to  achieve  the  highest 
accuracy  possible  in  the  comparisons  of  the  field  and  no  -  field  results. 

The  standard  deviation  of  average  samples  was  found  to  be  0.  025  pet. 
copper  which  gives  an  accuracy  in  the  difference  between  the  field  and  no  - 
field  segregation  curves  of  i  0.1  pet.  copper. 
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IV,  RESULTS 

A  photograph  of  a  typical  ingot  is  shown  in  Figure  8.  The 
normal  segregation  curves  for  the  4  inch  ingots  are  illustrated  in  Figures 
9  to  16.  Figure  17  shows  the  curves  for  the  2  inch  ingots.  The  field  and 
no  -  field  curves  are  shown  together  for  comparison.  Figures  10  and  ll 
show  the  type  of  segregation  that  occurs  if  there  is  entrapment  during  the 
initial  stages  of  freezing.  Figure  ll  is  typical  of  the  type  of  segregation 
obtained  for  a  relatively  fast  solidification  rate  (in  this  case  4  inches  per 
hour).  A.  macrograph  of  the  front  end  of  an  ingot  solidified  at  this  rate 
showed  dendritic  growth  and  entrapment  as  illustrated  in  Figure  18. 

Since  there  was  a  slight  amount  of  supercooling  and  entrapment 
at  the  beginning  of  solidification  in  each  ingot,  the  segregation  was  extrap¬ 
olated  to  the  front  of  the  ingot  in  order  to  obtain  comparative  effective 
distribution  coefficients.  The  method  of  extrapolation  balances  the  solute 
mass  across  the  enriched  (by  entrapment)  and  depleted  zones  of  the  ingot. 
The  values  of  the  effective  distribution  coefficients  corresponding  to  the 
different  rates  of  solidification  from  1/4  to  2  inches  per  hour  are  presented 
graphically  in  Figure  19.  The  K  values  are  calculated  using  the  relation¬ 
ship  Cj  =  KCc  for  the  entrapolated  concentration  at  the  ingot  end,  C0  was 
determined  by  balancing  the  solute  mass  across  the  whole  ingot. 


.  *  - 
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Figure  8.  Photograph  of  a  Typical  Unidirectionally  Solidified  Ingot. 


(full  size) 
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Normal  Segregation  in  the  A1  -  3  pet.  Cu  System  for  a  Solidification  Rate  of 
4  inches  per  hour  and  a  Magnetic  Field  of  34,  000  oersteds 
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Normal  Segregation  in  the  A1  -3.0  pet.  Cu  System  for  a  Solidification  Rate  of 
1  inch  per  hour  and  a  Magnetic  Field  of  34,  000  oersteds. 
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Figure  12.  Normal  Segregation  in  the  A1  -  3.0  pct.Cu  System  for  a  Solidification  Rate  of 
2  inches  per  hour  and  a  Magnetic  Field  of  34,  000  oersteds. 
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Normal  Segregation  in  the  A1  -  2.9  pet.  Cu  System  for  a  Solidification  Rate  of 
1  inch  per  hour  and  a  Magnetic  Field  of  34,  000  oersteds. 
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igure  14.  Normal  Segregation  in  the  A1  -2.8  pet.  Cu  System  for  a  Solidification  Rate  of 
l/2  an  inch  per  hour  and  a  Magnetic  Field  of  34,000  oersteds. 
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Normal  Segregation  in  the  A1  -  3.0  pet.  Cu  System  for  a  Solidification  Rate  of 
l/4  of  an  inch  per  hour  and  a  Magnetic  Field  of  34,  000  oersteds. 
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igure  17.  Normal  Segregation  in  a  two  inch  Ingot  of  A1  -  3  pet.  Cu  for  a  Solidification  Rate 
of  1  inch  per  hoar  and  a  Magnetic  Field  of  34,  000  oersteds. 
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Figure  18.  Macrograph  of  the  Front  of  an  Ingot  Showing 
Dendritic  Growth  and  Entrapment.  (x7) 
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V.  DISCUSSION 

Figure  19  shows  that  the  magnetic  field  significantly  increases 

the  amount  of  normal  segregation  in  aluminum  -  3  pet.  copper  alloys  at 

solidification  rates  exceeding  1/2  an  inch  per  hour.  The  field  solidified 

ingots  have  an  effective  distribution  coefficient  of  0.67  compared  to  0.76 

(61 

for  the  no  -  field  ingots.  Burton  et  al  have  shown  that  extrapolation  of 
this  curve  to  zero  rate  of  growth  yields  the  equilibrium  distribution  coeff¬ 
icient  since  £  approaches  zero  (cf.  Equation  4).  For  the  alloy  system 
studied  here,  the  equilibrium  distribution  coefficient  is  0.  17.  The  curves 
of  the  field  and  no  -  field  results  in  Figure  19  become  coincident  for 
solidification  rates  less  than  1/2  an  inch  per  hour,  and  appear  to  extrap¬ 
olate  to  the  equilibrium  value  of  0.  17  for  K0  at  zero  growth  rate.  The 
magnitude  of  the  effect  appears  to  level  off  at  a  freezing  rate  of  1  inch  per 
hour.  It  is  difficult  to  interpret  results  at  solidification  rates  exceeding 
2  inches  per  hour,  however,  because  the  mechanisms  of  normal  growth 
tend  to  break  down  and  dendritic  growth  becomes  predominant  (see  Figure 
10),  It  is  interesting  to  note  that  in  two  cases  where  a  substantial  amount 
of  supercooling  was  allowed  to  occur  in  the  first  liquid  to  solidify  (Figures 

9  and  10),  that  the  amount  of  inverse  segregation  was  increased  by  the 

(2) 

field.  This  agrees  with  the  results  obtained  by  Cahoon  '  '  for  inverse 
segregation  studies  of  field  -  cast  ingots.  The  field  does  not  affect  normal 
segregation  to  any  observable  extent  for  solidification  rates  less  than  1/2 


an  inch  per  hour. 


* 


. 
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It  remains  to  determine  the  process  by  which  the  magnetic 
field  affects  normal  segregation  for  freezing  rates  greater  than  1/2  an 
inch  per  hour.  According  to  the  relationship  of  Burton  et  al  (Equation  4), 
the  factors  which  determine  the  effective  distribution  coefficient  for  a 
given  freezing  rate  R,  are  the  liquid  diffusion  coefficient  D,  the  equilibrium 
distribution  coefficient  K0,  and  the  diffusion  layer  thickness  £ .  Thus,  it 
is  possible  to  decrease  K  by  adjusting  one  or  all  of  these  factors.  A 
combination  of  changes  is  suspected  since  the  effect  is  not  the  same  for 
all  freezing  rates.  The  hypothesis  put  forth  is  that  the  equilibrium  distrib¬ 
ution  coefficient  is  lowered  via  a  decrease  in  the  liquid  diffusion  rate  of 
the  solute  atoms.  The  decrease  in  K  also  implies  a  change  in  the  equil¬ 
ibrium  phase  relationship  of  the  alloy  system.  It  is  further  postulated 
that  this  change  is  a  shift  in  the  liquidus  line  of  the  phase  diagram  such  as 
that  shown  in  Figure  20  which  results  in  a  decrease  in  K0.  Cahoon  '  '  also 
found  that  it  was  possible  to  explain  his  inverse  segregation  results  by 
such  an  alteration. 

It  is  not  possible,  however,  to  determine  quantitatively  the 
extent  of  this  liquidus  shift  from  normal  segregation  results.  At  freezing 
rates  of  about  l  inch  per  hour  a  12  pet.  decrease  in  the  effective  distribution 
coefficient  was  observed  with  the  field,  ie.  from  0.  76  to  0.  67,  but  it  does 
not  necessarily  follow  that  an  equal  change  occurs  in  the  equilibrium  dis¬ 
tribution  coefficient,  K0,  This  is  because  K  is  also  dependent  on  the 
diffusion  coefficient  and  if  the  latter  is  decreased  as  suggested,  a  larger 
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Normal  Liquidus 
"Field"  Liquidus 
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ure  20.  Partial  Phase  Diagram  of  the  Al(rich)  -  Cu  System 


Showing  the  Shift  of  the  Liquidus  Line. 
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decrease  in  Kc  is  required  for  the  observed  decrease  in  K.  The  change 

in  K0  corresponding  to  the  change  in  K  will  also  depend  on  the  effect  of  a 

(18) 

magnetic  field  on  convective  mixing.  Gupta'  '  has  shown  that  convective 
flow  transverse  to  a  magnetic  field  is  severely  inhibited.  Inhibition  of 
convective  mixing  tends  to  increase  the  effective  distribution  coefficient 
because  the  thickness,  £,  of  the  interfacial  solute  layer  is  increased  (cf. 
Equation  4).  This  indicates  that  the  magnetic  field  decreases  the  equil¬ 
ibrium  distribution  coefficient  even  more  than  suggested  above. 

There  is  no  observable  "field"  effect  at  very  slow  freezing 
(7) 

rates  because  Dean  et  al  have  shown  experimentally  that  the  diffusion 
coefficient  of  the  solute  atoms  has  an  ever  diminishing  effect  on  the  effective 
distribution  coefficient  as  the  solidification  rate  decreases.  Therefore,  at 
slow  rates  any  effect  of  the  magnetic  field  on  the  diffusion  coefficient  will 
not  be  manifested  by  a  change  in  K.  This  explains  why  the  field  does  not 
have  any  observable  effect  on  the  effective  distribution  coefficient  for 
freezing  rates  less  than  1/2  an  inch  per  hour.  In  the  next  section  a  theory 
is  presented  which  qualitatively  justifies  the  shift  of  the  liquidus  of  the 
equilibrium  phase  diagram  of  the  alloy  system. 

It  was  considered  necessary  to  establish  whether  the  magnetic 
properties  of  the  component  atoms  of  the  alloy  studied  had  any  effect  on 
solidification  under  the  experimental  conditions.  Since  aluminum  is  para¬ 
magnetic  and  copper  is  diamagnetic,  these  metals  will  experience  opposing 
forces  when  situated  in  a  magnetic  field  gradient.  Thu^  normal  segregation 


. 
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may  have  been  influenced  by  the  sharp  field  gradient  in  the  liquid  metal 
outside  the  pole  pieces  of  the  magnet.  For  the  2  inch  ingot  run  the  liquid 
metal  was  entirely  enclosed  between  the  pole  pieces  and  thus  not  under  the 
influence  of  any  significant  magnetic  field  gradient.  Figure  17  shows  that 
the  field  still  had  the  same  effect  on  the  segregation  as  it  did  with  longer 
ingots  solidified  at  the  same  rate.  It  may  be  concluded,  therefore,  that 
magnetic  field  gradients  had  no  appreciable  effect  on  normal  segregation. 

The  segregation  curves  for  comparative  ingots  extrapolate  to 
the  same  concentration  at  the  ingot  end  where  solidification  was  initiated. 
This  is  in  good  agreement  with  the  postulated  theory  because  the  field  will 
not  affect  the  solidification  process  until  a  diffusion  layer  has  been  estab¬ 
lished  ahead  of  the  solid  -  liquid  interface.  A  finite  amount  of  solid  must 
be  formed,  however,  before  a  solute  concentration  gradient  can  be  built 
up.  Therefore,  the  field  should  not  be  expected  to  have  any  influence  on 
normal  segregation  in  the  initial  portion  of  the  ingot. 

Another  interesting  feature  of  the  results  is  the  pattern  of  the 
segregation  curves  in  the  portion  of  the  ingots  which  solidified  last.  In 
this  region  the  field  curves  tend  to  level  off  and  fall  below  the  no  -  field 
curves.  This  may  be  attributed  to  the  following  factors: 

(i)  As  shown  by  Tiller  and  Rutter,  there  is  a  tendency  for  a  planar 
interface  to  break  down  with  an  increasing  solute  concentration  in  the 
liquid.  This  results  in  dendritic  growth,  and  the  formalisms  of  normal 
solidification  can  no  longer  be  employed  to  describe  the  freezing  process. 


/  '' 
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(ii)  The  proposed  diffusion  inhibition  in  the  magnetic  field  results  in  a 
higher  solute  concentration  ahead  of  the  solid  -  liquid  interface. 

The  combination  of  these  effects  causes  solidification  in  the  field  to 
proceed  by  columnar  dendritic  growth  tending  to  give  a  macroscopically 
uniform  solute  concentration  throughout  the  region. 


; 
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VI.  THEORY  OF  MAGNETIC  EFFECTS 


In  the  discussion  of  the  results  it  was  proposed  that  the  magnetic 


field  effect  observed  on  normal  segregation  occurs  via  diffusion  inhibition  , 
with  a  consequent  change  in  the  position  of  the  liquidus  curve.  We  now 
consider  the  various  possible  effects  of  a  magnetic  field  and  show  that  the 
hypothesis  of  diffusion  inhibition  and  equilibrium  phase  relationship  change 
is  a  reasonable  one. 

A.  Thermodynamic  Considerations 
The  application  of  a  magnetic  field  to  a  solidifying  alloy  may 
affect  the  constitutional  diagram  of  the  system  through  a  change  in  the  heat 
of  solidification.  The  energy  of  a  unit  mass  of  material  of  specific  magnetic 

5 '  M  Q\ 

susceptibility  placed  in  a  magnetic  field  of  strength  H,  is  given  by'  ’ 


(5) 


The  change  in  the  heat  of  solidification  of  a  metal  solidified  in  a  magnetic 
field  will  thus  be 


(6) 


where  the  superscripts  refer  to  the  liquid  and  the  solid  respectively. 

The  magnetic  energy  is  small  for  weakly  diamagnetic  and  paramagnetic 
substances,  not  exceeding  a  few  calories  per  gram.  Using  aluminum  as  an 


example 

10 


the  energy  change  upon  .solidification 


emu  and  cgs  units  are  used  unless  otherwise  specified. 
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in  a  34,  000  oersted  field  is  6  x  10  calories  per  gram.  This  is  negligible 
compared  to  the  heat  of  solidification  of  94  calories  per  gram,  and  will 
have  no  observable  effect  on  the  constitutional  diagram  of  the  alloy  system. 

Another  magnetic  phenomenon  which  may  affect  the  thermo¬ 
dynamics  of  an  alloy  system  is  magnetostriction.  Magnetostriction  is 
significant  only  in  ferromagnetic  materials  and  will  not  produce  a  signif¬ 
icant  bulk  strain  energy  in  paramagnetic  or  diamagnetic  metals  large 
enough  that  its  effects  will  be  noticeable  in  segregation  work.  For  example 
A  l/l  for  bismuth,  a  highly  diamagnetic  metal,  is  only  of  the  order  of  10 
in  a  magnetic  field  of  IcA  oersteds.  This  very  small  distortion  of  the 

lattice  results  in  a  correspondingly  negligible  energy  change. 

Finally  it  is  necessary  to  consider  the  effect  of  a  magnetic 

field  on  the  solid  solubility  limits  of  an  alloy  system.  Phase  changes  in 

alloys  occurring  at  certain  electron/atom  ratios  are  attributed  to  a  con- 

122) 

traction  of  the  Fermi  sphere  which  lowers  the  electron  energy.  '  This 
occurs  when  the  Fermi  sphere  contacts  the  Brillouin  zone  boundary  and 
part  of  the  crystal  structure  changes  to  accomodate  more  solute  atoms 
without  increasing  the  electronic  energy.  If  the  magnetic  field  was  to 
change  the  electronic  energy  of  the  alloy  system,  a  particular  phase  change 
could  occur,  thus  changing  the  solid  solubility  limits.  A  magnetic  field 
coalesces  into  discrete  degenerate  levels  the  quasi  -  continuous  energy 
levels  of  the  conduction  electrons  in  a  metal,  the  degree  of  degeneracy 
depending  on  the  field  strength.  As  a  result  there  is  an  increase  in  the 
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mean  kinetic  and  exchange  energies  of  the  electrons.  Dresselhaus  '  ' 

has  shown  however,  that  the  electron  energy  increases  appreciably  only 
when  the  magnetic  energy  i^u)ce,  equals  the  Fermi  energy,  where  K  is 
Plank's  constant  andtbc^is  the  cyclotron  frequency  of  the  electron.  The 
cyclotron  frequency  is  given  by 


e_H 

m 


(7) 


where  e  is  the  charge  of  an  electron  and  m  is  its  effective  mass.  Since 
the  Fermi  level  for  metals  is  1  to  10  electron  volts,  a  field  strength  of 

Q 

about  10  oersteds  would  be  required  to  effect  any  appreciable  change  in 
the  Fermi  energy.  Thus,  this  mechanism  will  not  account  for  the  observed 
segregation  results,  and  we  must  look  elsewhere  for  an  explanation  of  the 
magnetic  effects  on  solidification. 


B.  Plasma  Concept  of  a  Metal 

In  its  original  and  narrower  sense  the  term  plasma  was  applied 

to  an  ionized  gas  having  a  nearly  equal  number  of  positive  and  negative 

charges  so  that  the  state  of  the  whole  is  electrically  neutral.  The  term 

plasma  is  now  generally  taken  to  mean  "the  portion  of  a  material  body 

which  is  much  larger  than  a  Debye  length  and  throughout  which  charged 

(24) 

particles  move".  '  The  Debye  length  is  a  shielding  distance  beyond 
which  the  electric  field  of  a  charged  particle  is  shielded  by  particles  of 
opposite  sign.  A  solid  or  liquid  metal,  therefore,  may  be  considered  to 
have  plasma  -  like  properties  by  virtue  of  the  presence  of  the  conduction 
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electrons.  Since  the  translational  freedom  of  the  ion  cores  in  the  liquid 
is  much  greater  than  that  of  the  ions  in  the  solid,  the  liquid  is  nearer  to 
the  original  concept  of  a  plasma  than  is  the  solid. 

The  study  of  the  interaction  of  magnetic  fields  with  plasmas 
or  moving  conducting  fluids,  if  the  continuum  property  is  present,  is 
known  as  magnetohydrodynamics.  It  is  necessary  to  examine  the  magnet¬ 
ohydrodynamic  forces  in  a  conducting  fluid  and  their  possible  effects  on 
the  transport  processes  of  solidification.  In  a  liquid  metal  under  the 
influence  of  a  magnetic  field  the  ensemble  of  electrons  and  ions  can  be 
viewed  as  charges  which,  move  across  the  magnetic  lines  of  force,  and 
they  thus  experience  electromagnetic  forces  that  accelerate  them  in  a 
direction  perpendicular  to  both  the  field  lines  and  their  instantaneous 
direction  of  velocity.  The  magnetohydrodynamic  equations  of  motion  for 
a  plasma  must  contain  quantities  such  as  pressure,  density,  viscosity, 
and  constraints  imposed  by  the  magnetic  field  variables.  Thus,  Maxwell's 
equations  must  be  coupled  with  the  hydrodynamic  equations  to  define  the 
motion  to  the  bulk  plasma. 

(25) 

The  magnetohydrodynamic  equation  of  motion  is  given  by 


p  T/  -f 

where: 


pfhJ-  V)TT  =  p|  Hr  T/  -h  M  j  A  H  , 


(8) 


v  =  macroscopic  velocity  of  the  fluid 


p  =  density  of  the  fluid 


1 
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p  z  pressure 

^  -  kinematic  viscosity  of  the  fluid 
yM  -  permeability  of  the  plasma 
ft  -  gravitational  acceleration 
H  =  magnetic  field  strength 
^  -  macroscopic  current  density 

Along  with  Equation  (8)  which  is  derived  from  the  Boltzmann  transport 

/  T  L  \ 

equation,  '  Maxwell's  equations  relating  the  electrical  and  magnetic 
fields,  E  and  H,  are  required  for  the  complete  description  of  the  system. 
However,  for  the  case  of  finite  resistivity  and  constant  magnetic  fields, 
the  system  may  be  described  by  Equation  (8)  and  the  equation  for  the 
macroscopic  current 


where  <T  is  the  electrical  conductivity  of  the  fluid.  The  electromagnetic 
force  x  H  in  Equation  (8)  is  often  referred  to  as  the  ponderomotive 
force.  It  tends  to  impede  motion  across  the  lines  of  force.  If  ^  is 
replaced  by  Equation  (9)  this  force  becomes 


- yu  Hx Vi ) 


(10) 


where  Ex  and  vx  are  the  components  of  E  and  v  transverse  to  the  magnetic 
field.  The  quantity ^(THv^  is  the  induction  drag  tending  to  destroy  motion 
across  the  lines  of  force,  and  is  regarded  as  a  magnetic  viscosity.  This 
induction  drag  will  inhibit  diffusion  in  the  liquid  only  to  the  extent  that 
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there  is  bulk  motion  of  the  fluid  associated  with  diffusion.  Diffusion  in 
liquid  metal  systems  is  relatively  slow  and  density  changes  are  small 
however,  so  the  effect  of  the  magnetic  viscous  force  on  solidification 
will  be  negligible. 

Although  the  macroscopic  forces  may  be  negligible,  the  magnetic 
forces  on  the  individual  particles  remain  and  will  affect  their  relative 
motion.  A  charged  particle  describes  between  collisions  a  spiral  path 
around  a  magnetic  line  of  force.  Equation  (8)  is  only  valid  when  the 
gyration  radius  of  a  particle  is  much  larger  than  its  mean  free  path.  If 
the  field  strength  is  such  that  the  gyration  radius  is  of  the  order  of  a 
mean  free  path,  the  system  must  then  be  described  by  the  fundamental 
kinetic  equations  of  plasma  dynamics. 

A.  pressure  gradient,  Vp,  in  a  plasma  under  the  influence  of  a 

magnetic  field  is  responsible  for  various  drift  motions  of  the  plasma 

perpendicular  to  the  field.  Neglecting  particle  collisions,  the  velocity  of 

this  drift  perpendicular  to  H  and  ^p  is  given  by'  ; 

-> 


(ID 


where  Vp.  is  the  pressure  gradient  of  the  ion  cores  and  ne  is  the  electron 
density.  If  collisions  between  electrons  and  ions  are  considered,  there 

^  /  7  Q  \ 

will  be  an  additional  contribution  to  ~vL  given  by'  ' 


- 

H2- 
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where  Vp  is  the  total  pressure  gradient  and  Oi  is  the  electrical  conduct¬ 
ivity  transverse  to  the  magnetic  field.  Thus,  for  liquids  of  finite  conduct¬ 
ivity  there  are  normally  two  drift  components,  one  parallel  to  the  pressure 
gradient  and  one  perpendicular  to  it.  The  component  of  vx  parallel. to  Vp 
represents  a  net  outward  drift  of  the  plasma.  Collisions  between  like 
particles  do  not  to  a  first  approximation  contribute  to  any  alteration  of 
the  density  distribution  since  an  encounter  between  two  identical  particles 

shifts  the  two  guiding  centers  (orbital  centers  of  the  gyrating  particles) 

( 29) 

by  equal  and  opposite  amounts.  ' 

"Diffusion"  caused  by  these  drift  motions  is  flow  resulting  from 
the  motion  of  ions  and  electrons  together,  not  movement  of  one  species 
relative  to  another.  It  is  therefore  not  possible  to  make  a  comparison  of 
the  flow  component  along  Vp  in  Equation  (12)  with  normal  diffusion  velocities 
in  liquid  metals.  Equation  (12)  does,  however,  show  the  nature  of  diffusion 

inhibition.  The  extent  of  flow  or  "diffusion"  parallel  to  the  density  gradient 

2. 

is  inversely  proportional  to  H  ,  causing  flow  along  the  gradient  to  vanish 
for  very  high  field  strengths. 

A  more  useful  insight  into  the  diffusion  mechanism  can  be 
obtained  if  the  alloy  is  viewed  basically  as  a  two  fluid  plasma  of  ions  and 
electrons.  The  relative  effect  of  the  field  on  the  cyclic  motion  of  the 
electrons  to  its  effect  on  the  ions  is  proportional  to  the  inverse  ratio  of 
their  masses.  The  motion  of  the  electrons  only  need  be  considered,  for 
through  the  requirement  of  electrical  neutrality,  the  electrons  will  determine 
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the  motion  of  the  ions  (ambipolar  diffusion).  A  concentration  gradient 

at  the  solid  -  liquid  interface  results  in  a  free  electron  gradient  (providing 

the  valencies  of  the  component  metals  are  different),  which  in  a  magnetic 

field  is  assumed  to  behave  as  a  plasma  density  gradient.  If  the  electrons 

are  able  to  describe  at  least  a  turn  or  so  of  their  spiral  paths  between 

successive  collisions,  their  average  displacements  along  the  density 

(30) 

gradient  will  be  greatly  reduced  by  the  presence  of  the  field.  If  this 

situation  is  realized  in  a  liquid  metal  in  a  moderate  magnetic  field,  the 

immoblization  of  electron  motion  across  the  field  would  in  turn  determine 

the  ion  motion  because  of  the  basic  requirement  of  electrical  neutrality. 

( 31) 

Assuming  that  the  zone  structure  is  largely  retained  in  the  liquid  state, 
the  calculated  value  of  the  gyration  radius  for  an  electron  moving  at  the 
Fermi  velocity  in  a  field  of  10  oersteds  does  approach  the  normal  mean 
free  path  of  the  electron.  The  extent  to  which  the  charge  distribution 
within  a  plasma  can  deviate  from  the  equilibrium  distance  is  the  Debye 
distance,  h^,  given  by 


(13) 


where  n£  is  the  electron  density,  k  is  Boltzmann's  constant  and  T  is  the 
absolute  temperature.  In  metals  this  distance  is  only  the  order  of  an 
angstrom  unit  if  the  thermal  energy,  kT,  is  replaced  by  the  Fermi  energy. 
Therefore,  the  heavy  positive  ions  will  be  confined  by  the  restrained 

electrons,  and  diffusion  will  be  inhibited. 

A  useful  approximation  for  evaluating  hDby  which  order  of  magnitude 
agreement  is  obtained  with  more  refined  theories  (Thomas -Fe rmi 
approximation  of  the  charge  potential  for  example) 
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C.  Irreversible  Thermodynamics  and  the  Diffusion  Constraint 

In  order  to  determine  the  effect  of  diffusion  inhibition  on  the 

solidification  process,  the  irreversible  thermodynamics  of  the  process  must 

be  considered.  The  thermodynamic  principle  governing  an  irreversible 

process  is  that  the  rate  of  entropy  production  is  positive  definite  and  a 
(32) 

minimum'  The  minimum  rate  of  entropy  production  may  be  described 

by  a  variational  principle,  where  variation  is  made  with  respect  to  the 

(33) 

free  thermodynamic  forces.  The  rate  of  entropy  production  per  unit 

volume  is  represented  as  a  bilinear  function  of  the  fluxes  and  their 

— ► 

conjugate  forces  Xt* .  For  solidification  of  a  binary  alloy  these  fluxes  and 

forces  are  those  related  to  the  temperature  gradient,  the  free  energy 

difference  for  the  liquid  -  solid  reaction,  and  the  chemical  potential 

(34) 

gradient.  Kirkaldy  has  shown  that  the  variational  principle  generates 

Le  Chatelier's  principle  of  moderation  in  the  presence  of  constraints  for 

steady  state,  heterogenous  convective  systems  such  as  crystal  growth  in 

binary  alloys.  The  variational  principle  can  also  be  used  to  define  the 

(35) 

process  path  of  non  -  stationary  heterogenous  dissipative  systems. 

The  entropy  production  rate  is  still  a  minimum,  but  it  must  now  be 
consistent  with  the  external  dynamic  constraints  that  assure  a  non  -  zero 
transformation  rate. 

The  system  under  consideration  has  an  undefined  internal 
constraint  associated  with  the  inhibited  diffusion  process  which  is  coupled 
in  some  unknown  way  to  the  interfacial  reaction.  This  internal  constraint 


- 
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may  be  replaced  by  a  suitable  external  mathematical  one  which  allows 

conventional  phenomenological  formalism  and  so  permits  definition  of  the 

(35) 

system  in  terms  of  variational  principles.  The  actual  process  path 

the  system  will  chose  will  be  that  which  is  compatible  with  the  dynamic 
constraints  and  which  tends  to  minimize  the  rate  of  entropy  production. 

de  Groot  gives  for  the  entropy  production  rate,  C  ,  of  a  non¬ 


polarized  continuous  system  undergoing  a  thermal  conduction  -  diffusion  - 


viscous  flow  process  in  a  electromagnetic  field  the  following  expression 


.(36) 


<T  =  (-  J,-,jradT  ^  v)/T  >0  (14) 


Equation  (14)  consists  of  a  sum  of  products  of  fluxes  and  thermodynamic 

forces  containing  terms  arising  from  heat  conduction,  diffusion,  electric 

/  *  — >  ’  * 

conduction,  and  viscous  flow  respectively.  Ef-E  +  vXH)  is  the  electric 
field  relative  to  the  axis  moving  at  the  velocity  v  of  the  mass,  TT  is  the 
viscous  pressure  tensor,  and  u^  represents  the  chemical  potentials  of 
the  diffusing  substances.  The  entropy  production  rate  differs  from  that 
without  an  electromagnetic  field  only  by  the  inclusion  of  the  electric 
conduction  term,  i*E.  This  is  a  Joule  heat  dissipation  term  arising  from 

Hall  currents.  The  effect  of  the  field  is  found  in  the  diffusion  and  the 

— > 

viscous  flow  term  since  it  is  assumed  that  JK  and  the  transverse  component 

— >  — > 

of  v  are  reduced.  The  rate  of  entropy  production  increase  i*E/T  should 
be  equal  to  the  decrease  in  the  diffusion  entropy  factor.  If  the  solidification 
time  for  the  system  is  unaltered  by  the  field,  the  diffusion  and  viscous  flow 
terms  will  decrease  the  total  entropy  production  in  the  system  because 
the  Joule  heat  will  be  dissipated.  Heat  conduction  should  play  the  dominant 
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role  in  energy  dissipation  for  solidification  in  a  magnetic  field  because 
diffusive  flow  is  inhibited  whereas  transverse  heat  flow  is  affected  consid¬ 
erably  less  than  mass  transfer  itself.  Non  -  convective  heat  flow  at  liquid 

metal  temperatures  involves  momentum  transfer  largely  via  vibrating 
(37) 

ions,  so  its  component  transverse  to  the  field  will  not  be  affected  to 

the  same  extent  as  would  mass  transfer.  The  decreased  diffusion  and 
viscous  flow  entropy  production  will  manifest  itself  as  a  chemical  (segre¬ 
gation)  phenomenon. 

D.  Consequences  of  Diffusion  Inhibition 
Having  proposed  that  the  magnetic  field  inhibits  diffusion,  the 
effect  of  this  phenomenon  on  the  concentration  gradient  ahead  of  an  advancing 
solid  -  liquid  interface  must  be  considered.  The  concentration  gradient 
ahead  of  an  advancing  interface  is  shown  schematically  in  Figure  21.  If  the 
magnetic  field  inhibits  diffusion  of  the  solute  atoms,  the  mass  transfer  will 
tend  from  one  of  diffusion  to  one  of  forced  convective  -  viscous  flow.  Thus, 
a  "pile  -  up"  of  solute  atoms  at  the  interface  will  result  and  the  solute 
concentration  will  change  from  the  normal  form,  ,  to  CF  shown  in 
Figure  21.  If  the  equilibrium  solidification  temperature  is  plotted  for  this 
"field"  concentration  curve  using  the  unchanged  phase  diagram,  a  curve 
similar  to  that  illustrated  by  T^  is  obtained.  An  equilibrium  solidification 
temperature  gradient  such  as  T^  is,  however,  incompatible  with  the  heat 
flow  requirements  necessitated  by  the  temperature  boundary  conditions. 
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-  Normal  concentration  gradient  in  the  liquid 
Cp  -  Proposed  concentration  gradient  in  the  liquid 
T*-  Solidification  temperature  for  Cp  using  the  unchanged 
constitutional  diagram 

Tp  and  TF  -  Solidification  temperatures  for  Cpobtained  by  raising 
liquidus  line  of  the  constitutional  diagram 

Figure  21.  Concentration  Gradients  and  Temperature  Distributions 
at  the  Solid  -  Liquid  Interface 
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The  temperature  distribution  at  the  sink  must  have  the  form  ci  T/bx  <0, 
as  this  is  the  only  solution  possible  for  the  heat  flow  equation.  If  the  liquidus 
line  takes  the  form  of  the  dashed  line  in  Figure  20,  the  constitutional  require¬ 
ments  coupling  the  thermal  and  diffusion  fields  result  in  a  temperature 
gradient  of  the  form  TF  or  ever  Tp  ,  ie.  one  that  is  compatible  with  the  heat 
flow  requirements  of  the  system.  Therefore,  in  order  for  solidification 
to  occur  under  the  imposed  diffusion  constraint,  the  constitutional  diagram 
must  be  altered  as  proposed  (see  Figure  20). 


E.  Energy  Considerations 

It  remains  to  reconcile  the  proposed  change  in  the  constitutional 
diagram,  which  we  suggest  is  effected  via  diffusion  inhibition,  with  the  local 
equilibrium  principle.  For  equilibrium  considerations  the  total  or  cohesive 
energy  of  the  system  must  be  examined.  The  cohesive  energy  of  a  lattice 
can  be  represented  through  the  Schrodinger  wave  equation,  using  approx¬ 
imation  methods  for  its  solution.  The  cellular  approximation  of  Wigner 
and  Seitz  has  been  applied  with  some  success  to  bi-  and  multi-valent  metals. 

The  cohesive  energy  of  a  lattice,  E-,.  ,  referred  to  its  sublimated  state  is 

,  (38) 

given  by 


Eo  + 


3_  ^ 


(15) 


where  EQ  is  the  zero  point  energy  of  the  electron,  Z  is  the  valence  of  the 
metal,  and  r0  is  the  radius  of  the  unit  cell  in  the  sphere  approximation.  The 
second  term  is  the  average  Fermi  energy  of  the  electrons  and  be*/ r0  accounts 


' 
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for  the  electron  interaction  energy  per  electron  pair  in  any  one  atomic 

i/3 

cell.  Since  rD  is  proportional  to  V0  ,  the  atomic  cell  volume,  Equation  (15) 
emphasises  the  relationship  between  the  cohesive  energy  and  the  specific 
volume  of  a  metal. 

In  a  binary  alloy  system  it  is  usually  assumed  that  the  free 

electrons  exist  in  two  sets  of  Fermi  energy  levels  associated  with  the 

(39) 

potential  fields  of  the  constituent  ions.  The  two  sets  of  energy  levels 

are  superposed  to  give  a  single  energy  band  of  non  -  uniform  density 
(providing  the  energy  levels  are  different).  There  is  a  redistribution  of 
electrons  between  bands  to  lower  the  total  energy  by  an  amount  approx¬ 
imately  equal  to  the  bond  energy  of  the  system.  This  phenomenon  of  Fermi 
energy  equilibration  also  occurs  when  two  different  systems  are  brought 
into  mutual  contact.  Varley^^  has  shown  that  the  radius  of  an  atomic 

cell  in  alloys  is  approximately  given  by 


where  n  (n<0)  is  the  charge  transfer.  Therefore,  redistribution  of  charge 
will  result  in  an  atomic  volume  change  which  will  in  turn  affect  the  total 
energy  of  the  system. 

The  proposed  change  in  the  constitutional  diagram  per  se 
represents  a  departure  from  local  equilibrium,  and  therefore  an  increase 
in  energy.  It  is  suggested  that  this  energy  increase  is  counterbalanced  by 
a  change  in  the  specific  volume  of  the  liquid  alloy  which  is  due  to  transfer 
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of  electrons  across  the  solid  -  liquid  interface  in  accordance  with  the 

principle  that  the  Fermi  levels  of  the  two  contiguous  phases  in  equilibrium 

must  be  coincident.  The  specific  volume  change  would  be  almost  entirely 

in  the  liquid  since  the  modulus  of  elasticity  of  the  liquid  is  about  1/3  that 

of  the  solid.  Any  volume  change  in  the  solid  would  increase  the  overall 

strain  energy  of  the  system.  Thus,  local  equilibrium  can  be  restored  at 

the  interface  through  the  process  of  electron  transfer  and  specific  volume 

increase  or  decrease  in  the  liquid.  This  volume  change  proposal  has  been 

(2) 

experimentally  supported  by  Cahoon.  '  He  found  that  it  was  possible  to 
account  for  an  increase  in  inverse  segregation  of  Al  -  Cu  alloys  in  a 
magnetic  field  only  by  altering  the  constitutional  -  specific  volume  relation¬ 
ship  of  the  system. 
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VII.  SUMMARY  AND  CONCLUSIONS 

The  present  investigation  of  the  effect  of  a  magnetic  field  on 
solidification  mechanisms,  principally  that  of  normal  segregation  and  the 
transport  processes  involved  therein,  leads  to  the  following  conclusions: 

(1)  Solidification  in  a  magnetic  field  increases  the  normal  segregation  in 
aluminum  -  3  pet.  copper  alloys  for  certain  freezing  rates.  The  amount 
of  segregation  for  rates  from  3/4  to  2  inches  per  hour  is  increased  by 

12  pet.  (the  effective  distribution  coefficient  is  reduced  from  0.76  to  0.67) 
and  below  rates  of  1/2  an  inch  per  hour  the  field  has  no  observable  effect. 
The  effective  distribution  coefficients  extrapolate  to  the  equilibrium  value 
at  zero  growth  rate. 

(2)  A.  theory  which  considers  the  plasma  -  like  properties  of  liquid  metals, 
and  which  is  based  on  the  effect  of  a  magnetic  field  on  irreversible  transport 
processes  is  developed  to  explain  the  results. 

(3)  The  theory  shows  that  diffusion  inhibition  caused  by  the  magnetic  field 
has  the  effect  of  raising  the  liquidus  curve  of  the  alloy  system.  This  is 
shown  to  be  compatible  with  heat  flow  requirements  for  the  solidifying 
system  in  the  presence  of  the  thermal  and  diffusion  constraint.  Local 
equilibrium  is  restored  at  the  interface  through  the  process  of  electron 
transfer  across  the  interface  and  the  accompanying  specific  volume  increase 
or  decrease  (in  this  case  an  increase)  in  the  liquid  alloy. 

(4)  Normal  segregation  results  do  not  give  quantitatively  the  extent  of  the 
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liquidus  shift  because  the  experimentally  determined  effective  distribution 
coefficients  are  not  directly  related  to  the  equilibrium  distribution  coeff-: 
icient.  The  effect  of  the  magnetic  field  inhibiting  convective  flow  and 
diffusion  will  tend  to  increase  K.  However,  a  decrease  of  12  pet.  was 
observed  in  the  effective  distribution  coefficient  for  solidification  in  the 
field.  This  indicates  that  the  magnitude  of  the  change  in  the  equilibrium 
distribution  coefficient  is  greater  than  the  decrease  of  the  effective  distrib¬ 
ution  coefficient.  Thus,  it  may  be  concluded  that  K0  for  Al  -  3  pet.  Cu 
alloys  is  decreased  somewhat  more  than  12  pet.  by  a  magnetic  field  of 
34,000  oersteds. 
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